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Abstract. Our goal is to assess Gaia's performance on the period re- 
0^ covery of short period (p < 2 hours) and small amplitude variability. 

. To reach this goal first we collected the properties of variable stars 

that fit the requirements described above. Then we built a database 
of synthetic light-curves with short period and low amplitude variabil- 
I i ity with time sampling that follows the Gaia nominal scanning law and 

with noise level corresponding to the expected photometric precision of 
Gaia. Finally we performed period search on the synthetic light-curves 
to obtain period recovery statistics. This work extends our previous 
period recovery studies to short period variable stars which have non- 
stationary Fourier spectra. 

t-H ■ 

> 

^ ■ 1 Introduction 

00 

^-j- ■ Variability is a common phenomena of stars. In the last decades the large pho- 
tometric surveys observed different fraction of variability, probably because these 
^ I . surveys had different photometric precision, time resolution and time sampling. 
C^, Gaia expected to detect 5-15 percent of variability (Eyer & Cuypers 2000) This 
t-H . means that around ~T00 million variable star light curves will have to be pro- 
cessed by the Gaia DPAC CU7. To have a glimpse about how Gaia will deal with 
• *h . variable see Eyer et al.'m these proceedings. 

Studying variable stars contribute to many field of astrophysics, just to mention 
• a few: With the methods of asteroseismology one can look inside the stars accessing 
information of the stellar interior which let one to probe stellar evolution and 
stellar pulsation theories. With the famous period-luminosity relations of classical 
Cepheids and RR Lyrae stars one can determine the distance scales inside the 
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Milky Way and to nearby galaxies. Studying compact pulsators like ZZ Ceti stars 
one can explore the physics of degenerate matter a regime which is not or hardly 
accessible in the laboratories. 

Our study is focusing on stars which shows variability in timescale shorter than 
2 hours, hereafter short period variables. We present their types and properties 
in Sect. [5] . These kind of variable stars, once detected, are advantageous from 
observational point of view as in a few hours long observation run many variation 
cycle is possible to be followed-up. In the scope of Gaia satellite, we are interested 
to know that what kind of information can be obtained from the observed Gaia 
light-curves of short period variables. For example detection of short period vari- 
ability is one of the questions we are interested in, the possibility to recover their 
periods is another one. In this proceedings we are exploring the latter case, by 
analysing the Fourier spectra of simulated Gaia light-curves of short period vari- 
able stars. We choose the type ZZ Ceti as basis for our analysis, because observed 
ZZ Ceti star light-curves show a few level of complexity which represents mostly 
the of observed behaviours of short period variable stars. We are discussing such 
observed complexities together with the types and properties in Sect. [SJ In the 
Sect. [3] we describe the mission properties of Gaia which are relevant to our study, 
in the Sect. [4] we describing the simulation and in the Sect. [5] we presenting the 
period recovery results on simulated data. 

2 Types and properties of short period variable stars 

Many type of variable stars show variability in time scales shorter than 2 hours as 
shown in the Table [T] Most of the types are non-radial pulsators so their variabil- 
ity amplitudes are in the millimagnitude level and they often shows multiperiodic 
variability. In some of the short period variable star types nonlincaritics in the 
light-curves are commonly observed. For example in the case of ZZ Ceti stars two 
different attempts were made to theoretically explain and model such nonlineari- 
ties. One of them is a model developed by Brickhill (1990)) then later presented, 
explained and used by others fWu fTMTl Goldreich & Wu lT990l Wu [20CTH Ising & 
Koester f200H Montogmery 2005) which includes interaction between the pulsation 
and convection zone, thus relates the origin of nonlinearities to the propagation 
of the flux trough the convective zone. According to the other model by Brassard 
et aL (|1995p the deviations from linearity in the emergent flux are caused by the 
radiation transfer process. These nonlinearities becomes larger if the amplitude 
of the temperature perturbation increases. This process is neglected in the theory 
of Brickhill, however models based on their calculation can be fitted well to the 
observed mono- and multiperiodic ZZ Ceti stars (Montgomery 2005; Montgomery 

mm . 

Further complexities are also observed in the ligth-curves of ZZ Ceti stars. 
For example Handler et al |2"DD ff reported amplitude and frequency variation of the 
modes of the ZZ Ceti star EC 14012-1446. The cause of the amplitude changes 
on timescales from weeks to years of white dwarf variables is not well understood 
(see Winget & Kepler EH]). 
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Table 1. Types and properties of short period variables 



Type 


Periods [minutes] 


Amplitudes [mag] 


/3 Cep stars 


96 - 480 


< 0.1 


S Scuti stars 


28 - 480 


0.003 - 0.9 


roAp stars 


6 - 21 


< 0.01 


EC14026 stars 


1.3 - 8.3 


< 0.03 


Betsy stars (PG1716) 


33 - 150 


< 0.03 


ZZ Ccti stars (DAV) 


0.5 - 25 


0.001 - 0.3 


V777 Her stars (DBV) 


2 - 16 


0.001 - 0.2 


GW Vir stars (DOV + PNNVs) 


5 - 85 


0.001-0.2 


Brown Dwarf pulsators 


60- 120 




eclipsing white dwarfs 


> 6 


< 0.75 



3 Relevant mission properties 

Gaia will observe about one billion stars with unprecedented astrometric and pho- 
tometric precision (see Fig. []} down to magnitude G 20 mag. Over its 5 year 
mission, the two field of view of Gaia will systematically scan all the sky according 
to the Gaia nominal scanning law and the observed sources will cross the CCDs 
of the common focal plane about 70 times on average. The detector in the focal 
plane consist 9 CCD in the astrometric field and each of them integrates charges 
for 4.4 seconds as the star crossing the CCD in TDI (Time Delay Integration) 
mode. This means 9 observation during one cross, and combining all the focal 
plane crosses a star on average will have about 630 CCD measurements in the G 
band. This is called G band per-CCD photometry, the averaged value of the 9 
CCD measurements called Field of View (FoV) transit photometry, the averaged 
value over all the mission called end of mission photometry. 

The time gap between two consecutive field of view transit is defined by the 
satellite design and by the Gaia nominal scanning law. This semi-regular time 
sampling defines the spectral window (SW) of a star on a given position on the 
sky. This spectral window usually has high values (SW(f) ~ 0.6 — 0.9) at higher 
frequencies (/ >= 12 [c/d]) but sometimes a higher peak (SW(f) ~ 0.6 )appears 
around 4 [c/d] which corresponds to the rotation of the satellite. 

The variability amplitudes of short period variables are mostly at millimagni- 
tude level. As comparison the expected photometric precisions of Gaia is also at 
the millimagnitude level in wide range of G magnitudes as it is shown in Fig. [1] 

4 Light-Curve simulations 

In order to simulate G band Gaia time series of short period variable stars, first we 
are using a light-curve model to generate continuous light-curves, then we sample 
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Gaia G magnitude [mag] 

Fig. 1. Expected photometric precisions of Gaia: per-CCD, field transit, and in the end 
of mission. Figure by M.Varadi using Gaia DPAC CU5 Photpipe. 

the continuous model light-curve with the nominal scanning law of Gaia, in the end 
we add Gaussian noise which corresponds to the expected photometric precision of 
Gaia (see Fig.Q}. We use the following models according to the required complexity 
level: mono-periodic, multiperiodic, multiperiodic with nonlinearities (ZZ Ceti), 
multiperiodic with non stationary spectra. The inputs of the models are periods 
and amplitudes which corresponds to their observed values in real ZZ Ceti stars. 

5 Analysis of period recovery statistics 

We derive the period recovery statistics from a sample of 200 light-curves. Each of 
them was generated with equal number of measurements, with different time sam- 
pling (these corresponds to stars with different sky positions) and with the same 
periods and amplitudes as inputs for the models. Then we compute the Deeming 
periodogram (see Deeming 1975) for each light-curve until a chosen maximum fre- 
quency value. After this we select the maximum peaks frequency and we compare 
with one of the input periods. If there is a match within the width of the peak 
then we count the period as successfully recovered. After pre-whitening we repeat 
the process if we are interested in the recovery of secondary frequencies. 

5.1 Mono-periodic model 

Eyer and Mignard 2005] made an extensive research on the correct detection rate 
of mono-periodic signals of Gaia Field of View transit time series for a wide range 
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Fig. 2. Period recovery percentages of mono-periodic signals in the short period regime. 
Solid lines corresponds to per-CCD photometry, Dashed lines corresponds to Field of 
View transit photometry. The signal to noise ratio is equal to 2.5 on the per-CCD 
measurements. 



of periods. They concluded that periods of regular variable star can be recovered 
even from signals with low S/N ratio and that the period recovery depends mainly 
on the ecliptic latitude. In the Fig. [2] we present a result which extends their 
work to the short period regime and compares the period recovery from per-CCD 
photometry with period recovery from Field of View transit photometry. The 
figure shows that for periods smaller than 5 minutes the period recovery on the 
per-CCD photometry is significantly higher than the period recovery on Field of 
View transit photometry. This is because the duration of the transit starts to be 
comparable with the period of the signal and in the case of per-CCD photometry 
this transit duration which is now significant part of the period is well sampled 
while it is averaged to a single value in the case of Field of View photometry. 

5.2 Multiperiodic models 

The multiperiodic models consist three cases: The first model consists a sum 
of sinusoidal signals with periods and amplitudes corresponding to the observed 
periods and amplitudes of a ZZ ceti star. The second is based on the first but taking 
into account the nonlinearities caused by the convection zone as we described in 
Sect. [SJ The nonlinear light-curve simulator code is the same what we described 
in Varadi et al 20051 The third model is composed from two models of the first 
kind, let us call them A and B. The model A is computed on the first half of the 
dataset (the first 2.5 year per-CCD observations) with a given input periods and 
amplitudes and the model B is computed on the second part of the dataset with 
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different input periods and amplitudes. The input periods and amplitudes for the 
models A and B were taken from two different observational epoch of the ZZ Ceti 
star analysed by Handler et o/[JT)08, This way we simulated Gaia time series of a 
ZZ Ceti star with non-stationary spectra. 

We already published our results on the period recovery percentages on for the 
cases with first and second model (see Varadi et aZ.2009 ). Now we extend it with 
period recovery percentages for the model with the highest level of complexity. We 
found that for a G = 15 magnitude star about 73.5 percent of the cases the the 
correct period can be recovered if we use model of the first kind. This percentages 
drops only 1.5 percent if we use the second kind of model but it drops down to 
17.5 percent if we use simple model with non-stationary spectra as described by 
the third kind of model. 

We conclude that the non-stationarity of the spectra which is a commonly 
observed feature of white dwarf variables can cause serious challenge to our period 
recovery attempts on white dwarf variables from Gaia time series. Fortunately 
subdwarf variables (the types EC 14026 stars and Betsy stars) are having more 
stable modes which keeps their spectra stationary many years. (Charpinet 2009 
private communications) . 
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